Determination of the number of if)' events at BESIII 
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The number of tp' events accumulated by the BESIII experiment from March 3 through April 14, 
2009, is determined by counting inclusive hadronic events. The result is 106.41 x (1.00±0.81%) x 10® . 
The error is systematic only; the statistical error is negligible. 



I. INTRODUCTION studies of V'' decays. The number of ip' events, N^> , is 

important in all ip' analyses, including studies both of 

In 2009, the world's largest tp' sample to date was col- 
lected at BESIII, allowing more extensive and precise 
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the direct decays of the ip' , as well as its daughters, Xcj, 
/ic, and rjc- The precision of N^i will directly affect the 
precision of all these measurements. 

In this paper, we determine N^i with -0' 
inclusive hadrons, whose branching ratio is known 
rather precisely, (97.85 ± 0.13)% Also, a large off- 
resonance continuum data sample at Ecm — 3.650 GeV 
with an integrated luminosity of 44 pb~^ was collected. 
These events are very similar to the continuum back- 
ground under the ifj' peak. Since the energy difference is 
very small, we can use the off-resonance data to estimate 
this background. 

BEPCII is a double-ring e+e^ collider designed to 
provide e+e^ interactions with a peak luminosity of 
10'^'^ cm~^s~^ at a beam current of 0.93 A. The cylin- 
drical core of the BESIII detector consists of a helium- 
based main drift chamber (MDC), a plastic scintillator 
time-of- flight system (TOF), and a CsI(Tl) electromag- 
netic calorimeter (EMC), which are all enclosed in a su- 
perconducting solenoidal magnet providing a 1.0 T mag- 
netic field. The solenoid is supported by an octagonal 
flux-return yoke with resistive plate counter muon iden- 
tifier modules interleaved with steel. The acceptance for 
charged particles and photons is 93% over An stereo an- 
gle, and the charged-particle momentum and photon en- 
ergy resolutions at 1 GeV are 0.5% and 2.5%, respec- 
tively. 

The BESIII detector is modeled with a Monte Carlo 
(MC) simulation based on GEANT4 [3, For the sim- 
ulation of inclusive ip' decays, we use the evtgen gen- 
erator [ij. Known ip' decay channels are generated ac- 
cording to branching ratios in the PDG [ij; the remain- 
ing unknown decays are generated by the lundcharm 
model 0. 



II. EVENT SELECTION 

There are many types of events in the data collected 
at the V'' energy point, including tp' hadrons and lep- 
ton pairs (e+e~, fJ,'^n~, and t+t~), radiative returns to 
the J/ip, and J/tp decays from the extended tail of the 
J/ip Breit-Wigner distribution. In addition, there are 
non-resonance (QED) processes, which make up the con- 
tinuum background, including e^e^ — > 7* — >■ hadrons, 
lepton pairs, and e~^e~ e'^e^ -|-X (X=hadrons, lepton 
pairs). Non-collision events include cosmic rays, beam- 
associated background, and electronic noise. The signal 
channel is the process ip' — )■ hadrons. The data collected 
at the off-resonance energy include all of the above except 
tjj' — > hadrons and lepton pairs. 

Event selection includes track level selection and event 
level selection. At the track level, good charged tracks 
are required to pass within 1 cm of the beam line in the 
plane perpendicular to the beam and within ±15 cm from 
the Interaction Point (IP) in the beam direction. Photon 
candidate showers reconstructed from the EMC barrel 
region (| cos9\ < 0.8) must have a minimum energy of 25 



MeV, while those in the end-caps (0.86 < | cos 6*1 < 0.92) 
must have at least 50 MeV. The showers in the angular 
range between the barrel and end-cap are poorly recon- 
structed and excluded from the analysis. Requirements 
on the EMC cluster timing are applied to suppress elec- 
tronic noise and energy deposits unrelated to the event. 

At the event level, at least one good charged track is 
required. If the number of good charged tracks is larger 
than 2, i.e. Ngood > 2, no additional selection is needed. 
If Ngood — 2, where the Bhabha and dimuon events are 
dominant backgrounds, the momentum of each track is 
required to be less than 1.7 GeV/c, and the opening an- 
gle between the two tracks is required to be less than 
176° to suppress these backgrounds. Figures [T] and [5] 
show scatter plots of the momentum of one track ver- 
sus the momentum of the second track for MC simulated 
Bhabha events and inclusive MC events with two charged 
tracks, respectively. Figures [3] and S] show the open- 
ing angle distributions of MC simulated Bhabha events 
and inclusive MC events with two charged tracks, re- 
spectively. In addition, EmsiUel Ecm > 0.4 is required 
to suppress low energy background (LEB), comprised 
mostly of e+e^ — ^ e+e^ + X and double ISR events 
(e+e~ -)■ IisrIisrX). Here, Eligible denotes the visible 
energy which is defined as the energy sum of all charged 
tracks (calculated with the track momentum and assum- 
ing a TT^ mass) and neutral showers, and Ecm denotes 
the center-of-mass energy. Figure [S] shows the Eyisme 
distribution for data and inclusive MC events with two 
charged tracks. The excess in data at low energy is from 
the LEB. 

If Ngood = 1, at least two additional photons are re- 
quired in an event. From all photon pair combinations, 
the combination whose invariant mass, M^^, is closest to 
the 7r° mass is selected, and IM^y-M^^ol < 0.015 GeV/c^ 
is required. EmsiUe/ Ecm > 0.4 is also required to sup- 
press the LEB. Figure [S] shows the M^^ distributions in 
the tt" mass region for data and MC simulation. Figure [7] 
shows the EmsiUe distribution for data and inclusive MC 
events. The excess in data at low energy is from LEB. 



The average Z-direction vertex for an event is defined 

as 

V7 = — 

^ ^ good 

where Vz is the distance along the beam direction of 
the point of closest approach of a track to the IP. Fig- 
ure [5] shows the Vz distribution for ip' data after the 
above selection. Events satisfying \Vz\ < 4.0 cm are 
taken as signal, while events in the sideband region 6.0 
cm< \Vz\ < 10.0 cm are taken as non-collision back- 
ground events. The number of observed hadronic events 
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FIG. 1. The distribution of P2 versus Pi from MC simulated 
Bhabha events. The horizontal and vertical lines show the 



selection requirements to remove Bhabha and e e 



events. 




FIG. 2. The distribution of P2 versus Pi from inclusive MC 
events with two charged tracks. The horizontal and vertical 
lines show the selection requirements to remove Bhabha and 
e^e~ — >■ events. 



{N ) is determined by 



l\jobs. 7W- Ty- 

— signal sideband- 



(1) 



Another method to determine the number of hadronic 
events (described below) is to fit the average Z-vertex 
with a double Gaussian to describe the signal and a poly- 
nomial to describe the non-collision events. 



III. BACKGROUND SUBTRACTION 

In principle, the number of QED events can be esti- 
mated from: 



7V«^^ = C-a-e, 



(2) 



where C is the luminosity, and a and e are the cross- 
section and efficiency, respectively, a is usually obtained 
from theoretical prediction, and e is determined from MC 
simulation. 

However in this analysis, we use the large sample of 
off-resonance data collected at 3.65 GeV to estimate the 
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FIG. 3. The distribution of angle between tracks for MC simu- 
lated Bhabha events. The arrow shows the angle requirement 
used to remove most Bhabha events. 
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FIG. 4. The distribution of angle between tracks for inclusive 
MC events with two charged tracks. The arrow shows the 
angle requirement used to remove most Bhabha events. 



continuum background. The events remaining, after im- 
posing the same selection criteria in the off-resonance 
data, also form a peak in the Vz distribution, as shown 
in Figure [S) The same signal and sideband regions are 
used as for the ■0' data to determine the collision and non 
collision events. With this method, the continuum back- 
ground subtraction is independent of MC simulation, and 
little systematic bias is introduced. 




FIG. 5. The Evisibie/ Ecm distribution for Ngood ~ 2 events. 
Dots with error bars are data; the histogram is MC simula- 
tion, normalized to E^isiUe/ Ecm > 0.4. 
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FIG. 6. The 77 invariant mass (M-yj) distribution in the tt" 
mass region for Ngood ~ 1 events. Dots with error bars are 
data; the histogram is MC simulation. 



FIG. 9. The average Z vertex (Vz) distribution of hadronic 
events in off-resonance data. The curves are a double Gaus- 
sian to describe the signal and a polynomial to describe the 
non-collision events. 
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FIG. 7. The EvisiUc/Ecm distribution for Ngood = 1 events. 
Dots with error bars are data; the histogram is MC simula- 
tion, normalized to Eyisibie/Ecm > 0.3. 



The contributions from radiative returns to J/ip and 
J/ip decays from the extended tail of the Breit-Wigner 
are very similar at the ip' peak and ofF-resonance energy 
due to the small energy difference. They are estimated to 
be 1.11 and 1.03 nb at the ip' peak and the ofF-resonance 
energy point, respectively, and according to MC simula- 
tion, the efficiencies for the known continuum processes 
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FIG. 8. The average Z vertex (Vz) distribution of hadronic 
events in tp' data. The curves are a double Gaussian to 
describe the signal and a polynomial to describe the non- 
collision events. 



at the two energy points are also similar. Therefore, 
the off-resonance data can be employed to subtract both 
the continuum QED and Jf-ip decay backgrounds using 
a scaling factor, /, determined from the integrated lumi- 
nosity multiplied by a factor of 1/s (s = i?cm) to account 
for the energy dependence of the cross-section: 

f = 4^-^2 - 3.677, (3) 
£3.65 3.686^ 

where, Cjp' and £3.65 are the integrated luminosities for 
-0' data and 3.65 GeV data, respectively. 

The luminosities at the two different energy points 
are determined from e+e^ — >■ 77 events using the same 
track and event level selection criteria. At the track 
level, no good charged tracks and at least two showers 
are required. The energy for the most energetic shower 
should be higher than 0.7 x Ebeam while the second 
most energetic shower should be larger than 0.4 x Ebeam, 
where Ebeam is the beam energy. At the event level, 
the two most energetic showers in the ip' rest frame 
should be back to back, and their phi angles must satisfy 
178° < 101-021 < 182.0°. The luminosity systematic er- 
rors nearly cancel in calculating the scaling factor due to 
small energy difference between these two energy points. 
The / factor can also be obtained using luminosities de- 
termined with Bhabha events. It is found to be 3.685. 

Also of concern is the LEB remaining in the 0' events 
after the EyisiUe/ Ecm requirement. In order to test if 
the continuum background subtraction is also valid for 
these events, candidate LEB events are selected by re- 
quiring Eyisibie/ Ecm < 0.35 whcrc there are few QED 
events expected. Figures [10] and [11] show the compari- 
son of Eyisihie/ Ecm between peak and off- resonance data 
for Ngood = 1 and Ngood = 2 events, respectively. The 
agreement between the two energy points is good for 
these events. The ratios of the numbers of peak and off- 
resonance events for Ngood — 1 and Ngood = 2 are 3.3752 
and 3.652, respectively. Compared with the scaling fac- 
tor obtained from luminosity normalization in Eq.([3]), 
a difference of about 10% is found for Ngood — 1 while 
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TABLE I. iV°''=' for peak and off-resonance data (xlO^), and 
the detection efficiency for inclusive tj)' decay events deter- 
mined with f 06 X f 0® ip' — >■ inclusive MC events. 





Ngood > 1 


Ngood > 2 


Ngood > 3 


Ngood > 4 


ip' data 


f06.928 


102. 79f 


81.158 


63.063 


off-resonance data 


2.f92 


1.98 


0.704 


0.433 


e(%) 


92.9f2 


89.860 


74.624 


58.188 



FIG. 10. Comparison of LEB events between peak and 
off-resonance data for Ngood = 1 events. Dots with error bars 
are ^p' data, and the histogram is off-resonance data. 
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FIG. 12. The cos^ distribution for charged tracks. Dots with 
error bars are data; the histogram is MC simulation. 



FIG. 11. Comparison of LEB events between i/)' peak and off- 
resonance for Ngood = 2 events. Dots with error bars denote 
i/i' data, and the histogram denotes off-resonance data. 



there is almost no difference for Ngood = 2 events. These 
differences will be taken as systematic errors. 



The small numbers of events from -0' ^ e+e 



and T~^T~ in data that pass our selection do not need 
to be explicitly subtracted since ip' — > lepton events are 
included in the inclusive MC and those passing the se- 
lection criteria wiW contribute to the MC determined ef- 
ficiency, so that their contribution cancels. 

Table H] shows the number of observed hadronic events 
for different multiplicity requirements for ip' and off- 
resonance data. Figures [T^ [T^ and [TUshow the cos 9, 
Evisibiei and charged-track multiplicity distributions af- 
ter subtracting background. 



IV. NUMERICAL RESULT 



The number of ip' events is determined from 



peak 



f . Afobs. 

J of f — resonance 



(4) 



where, N"^^^^. is the number of hadronic events observed 
at the Ip' peak from Eq. (H]), N°ff_^^^^^^^^^ is the number 
of hadronic events observed at the off-resonance energy 




FIG. 13. The visible energy distribution. Dots with error 
bars are data; the histogram is MC simulation. 
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FIG. 14. The charged-track multiplicity distribution. Dots 
with error bars are data; the histogram is MC simulation. 
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TABLE II. N^i (x 10^) for diflerent diarged-track multiplicity 
requirements. 



N^,, 106.414 106.279 105.289 105.643 



point, Ecm = 3.650 GeV, with the same selection criteria 
as those for peak data, and e is the selection efficiency 
obtained from the inclusive V'' MC sample, the branching 
fraction of ip' inclusive hadron is included in the 
efficiency. The relevant numbers are listed in Table |T] for 
different Ngood selection requirements. The factor / is 
the scaling factor which has been introduced in Eq. (2). 
With these numbers, we obtain the numerical result for 
Nji,/ listed in Table HIl for different choices of Ngood- We 
take the result for Ngood > 1 as the central value of our 
final result. 



V. SYSTEMATIC UNCERTAINTIES 

The systematic uncertainties include the uncertainties 
caused by tracking, the event start time (Tq), trigger ef- 
ficiency, background contamination, the selection of the 
signal and sideband regions, etc. 



A. Tracking 

Generally, the tracking efficiency for MC events is 
higher than that of data according to various studies . 
Assuming the average efficiency difference between data 
and MC is 1% per track, the effect can be measured by 
randomly tossing out 1% of MC simulated tracks. Only a 
difference of 0.03% on iV^/ is found for Ngood > 1 events 
with and without this tracking efficiency change; N^' is 
not sensitive to the tracking efficiency. 



B. Charged-track multiplicity 

Figure [U shows that the MC does not simulate the 
charged-track multiplicity very well. The error due to 
charged-track multiplicity simulation can be estimated 
by an unfolding method, which is described as follows. 
The generated true charged multiplicity in MC simula- 
tion is even, i.e., 0, 2, 4, 6, 8, • • • . The observed MC 
multiplicity distribution is obtained after simulation and 
event selection. For example, if the generated true mul- 
tiphcity is 4, the observed multiplicities are 0, 1, 2, 3, 
or 4 with different probabilities. Therefore, an efficiency 
matrix, , which describes the efficiency of an event gen- 
erated with j charged tracks to be reconstructed with i 
charged tracks, is obtained from MC simulation. The dis- 
tribution of the number of observed charged-track events 
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FIG. 15. The distribution of total momentum for Ngood = 2 
events. Dots with error bars are data; the histogram is MC 
simulation. 

in data, N°''^' , is known. The true multiplicity distribu- 
tion in data can be estimated from the observed multi- 
plicity distribution in data and the efficiency matrix by 
minimizing the x^. The is defined as 

10 

X ^ ATobs. ' 

1=1 i 

where the Nj (j = 0, 2, 4, 6, 8, 10) describe the true 
multiplicity distribution in data and are taken as floating 
parameters in the fit. The simulation is only done up 
to a true multiplicity of 10, since there are few events 
at high multiplicity. The total true number of events 
in data can be obtained by summing all fitted Nj] it is 
105.96 X 10^ which is lower than the nominal value by 
0.4%. We take this difference as the error due to the 
charged-track multiplicity distribution. 



C. Momentum and opening angle 

For Ngood = 2 events, momentum and opening an- 
gle requirements are used to remove the huge number 
of Bhabha events. When the momentum requirement 
is changed from P < 1.7 GeV/c to P < 1.55 GeV/c, 
the corresponding M°'"'- for peak and resonance data, 
as well as the efficiency change, but the change in 7V^/ 
is only 0.05%. When the angle requirement is changed 
from 9 < 176° to 61 < 160°, the change in N^, is 0.01%. 
Therefore, the total uncertainty due to momentum and 
opening angle requirements is 0.05%. Figures [T5] and [TBI 
show comparisons between data and MC simulations for 
momentum and opening angle distributions after back- 
ground subtraction, respectively. 



D. LEB background contamination 

N^i is insensitive to the visible energy require- 
ment. The difference between a tight requirement. 
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FIG. 16. The distribution of opening angle between tracks 
for Ngood = 2 events. Dots with error bars are data; the 
histogram is MC simulation. 



Evisibie/ Ecm > 0.45, and no requirement is only 0.1%. 
Conservatively, an error of 0.1% is assigned due to the 
background contamination. 



E. Determination of number of hadronic events 

Two methods are used to obtain N°'"'- . The first is to 
directly count the numbers of events in the signal and 
sideband regions; the second method is to fit the Vz dis- 
tribution with a double Gaussian for the signal and a 
polynomial for the background, as shown in Figs. [5] and IH) 
A difference of 0.28% is found between these two meth- 
ods which is taken as the error due to the uncertainty 
from the iV°^*- determination. 



F. Vertex limit 

When Vr < 1 cm is changed to 1^ < 2 cm, N^/ changes 
by 0.35%, while if \Vz\ < 10 cm is changed to \Vz\ < 15 
cm, there is almost no change. Therefore, the difference 
of 0.35% is taken as the error from the vertex require- 
ment. 



G. Scaling factor 

The scaling factor can be obtained for two different 
QED processes, e+e^ — )• 77 and e~^e~ e+e~. The 
corresponding results are 3.677 and 3.685. The dif- 
ference on N^/ due to the / factor can be calculated 
by A/ ■ 7V^''^^;^^>i (3.650 GeY)/N^> = (3.685 - 3.677) • 
3. 1808/106. 3T= 0.023%. The slight difference indicates 
the uncertainty caused by the normalization factor is neg- 
ligibly small. 



H. Choice of sideband region 

We take \Vz\ < 4.0 cm as the signal region and 
6 < \Vz\ < 10 cm as the sideband region. A differ- 
ence of 0.45% in is found by shifting the sideband 
region outward by 1.0 cm, which is about la of the Vz 
resolution, i.e., the sideband region is changed from 6 
cm< \Vz\ < 10 cm to 7cm < jV^I < H cm. We take this 
difference as the error due to the uncertainty caused by 
choice of the sideband region. 



I. tt" mass requirement 

This requirement is only used for Ngood = 1 events. 
AT^/ has a slight change of 0.11% when the mass win- 
dow requirement is changed from \Mj~^ — M^o| < 0.015 
GeV/c^ to \M^^ -M^o\< 0.025 GeV/c^; this difference 
is taken as the uncertainty due to tt*^ mass requirement. 



J. The cross section of e+e t~^t 

Since the off-resonance energy point is not very far 
from TT threshold, cr(e"'"e~ — >■ t'^t~) does not vary as 
1/s between the off-resonance energy and the peak, 
as other QED processes. The difference between the ob- 
served and the cross section assuming a 1/s dependence 
causes a change of 0.17% in N^i. This change is taken 
as a systematic error. 



K. X + J/V) 

The MC assumes B{tp' X + J/tp) « 57% from 
the PDG while the CLEO experiment determined a 
branching ratio of 62% Using CLEO's result, a new 
inclusive MC sample was generated. The corresponding 
efficiencies are 92.912%, 89.761%, 74.838% and 58.528% 
for Ngood > 1, 2, 3 and 4, respectively. Compared with 
numbers in Table [H the efficiency differences between 
these two MC samples are negligible. 



L. Event start time determination 

The Event Start Time (EST) algorithm is used to de- 
termine the common start time of the recorded tracks 
in an event. The efficiency of the EST determination 
affects the resolution of tracks from the tracking algo- 
rithm. These efficiencies for different charged tracks, 
e, /i, TT, K, and p, and photons are studied with dif- 
ferent control samples for both data and inclusive MC 
events, for example, J/ip — > 7r"*"7r~7r'', TT~^Tr~pp, and 
-0' — > TT^Tr~ J/i/j, J/4> — l^l^, etc. All comparisons in- 
dicate that the efficiencies of the EST determination are 
high for both track and event level (> 98%) selection, and 
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the difference between data and MC simulated events is 
quite small 0.1%). We take this difference as the un- 
certainty caused by the EST determination. 



M. Trigger efficiency 

The fraction of events with Nngood > 2 is about 97%. 
The trigger efficiency for these events is close to 100.0% 
according to a study of the trigger efficiency @. For 
Ngood = 1 events, an extra tt" is required, and the hadron 
trigger efficiency for this channel is 98.7% .9]. Since the 
fraction of Ngood — 1 events is only about 3%, the uncer- 
tainty caused by the trigger is negligible. 



N. The missing 0-prong hadronic events 

A detailed topology analysis is performed for Ngood = 
events in the inclusive MC sample. Most of these 
events come from known decay channels, such as ip' — > 
X + J/ip {X — 77, 7r°7r°, and tt+tt"), ^p' JXcj, and 
-0' — i> e^e^, fj,^ . The fraction of pure neutral events 
is less than 1.0%. For the known charged decay modes, 
the MC simulation works well according to many com- 
parisons between data and MC simulation in Section 3. 
To investigate the pure neutral channels, the same se- 
lection criteria at the track level are used. The criteria 
at the event level include Ngood = and N^ > 3. The 
latter requirement is used to suppress e+e" 77 and 
beam-associated background events. The same selection 
criteria are imposed on the off-resonance data. FiguresflTl 
and [TSlshow the distribution of total energy in the EMC 
for data and inclusive MC events. The peaking events 
correspond to the pure neutral candidates, and the num- 
ber of events is extracted by fitting. The difference in 
the number of fitted events between data and inclusive 
MC events is found to be 17%. Therefore, the uncer- 
tainty due to the pure neutral events should be less than 
17% X 1% = 0.17%, and this is taken as the systematic 
error on the missing 0-prong events. 



O. B{xl>' hadrons) 

The uncertainty of — > hadrons) is very small 
according to the PDG [1], 0.13%, which is taken as the 
error due to uncertainty of -0' decays to hadronic events. 



P. Total error 

Table IIIII lists all systematic errors. The total sys- 
tematic error is determined by the quadratic sum of all 
errors. 




Eemc (GeV) 

FIG. 17. The distribution of total energy in the EMC with 
Ngood = for data. The dot-dashed line denotes the signal 
shape of 1/;' — > neutral channel, the dashed line denotes 
the background shape from QED processes, and the shaded 
region is the background shape from xp' decay. 




FIG. 18. The distribution of total energy in the EMC with 
Ngood = for inclusive MC events. The dashed line denotes 
the signal shape of 1/'' — >■ neutral channel, and the shaded 
region is the background shape from ^p' decay. 



TABLE III. The systematic error (%) 



Source 


Error 


Background contamination 


0.10 


A'^" determination 


0.28 


Choice of sideband region 


0.45 


Vertex selection 


0.35 


Momentum and opening angle 


0.05 


Scaling factor (/) 


0.02 


0-prong events 


0.17 


Tracking 


0.03 


Charged-track multiplicity 


0.40 


CT(e^e~ r'^r") 


0.17 


B{^P' + J/xP) 


0.00 


tt" mass requirement 


0.11 


EST determination 


0.10 


Trigger efficiency 


negligible 


B{ip' — >■ hadron) 


0.13 


Total 


0.81 
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VI. SUMMARY 

The number of ip' events is determined using 

xjj' — >■ hadrons. The large off-resonance data sample at 
l^cm = 3.65 GeV is used to estimate the background un- 
der the il>' peak. The number of '0' events taken in 2009 
is measured to be (106.41 ± 0.86) x 10^, where the error 
is systematic only and the statistical error is negligible. 
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